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ABSTRACT 

Tidal friction in cxoplanct systems, driven by orbits that allow for durable nonzero eccentricities at 
short heliocentric periods, can generate internal heating far in excess of the conditions observed in our 
own solar system. Secular perturbations or a notional 2:1 resonance between a Hot Earth and Hot 
Jupiter can be used as a baseline to consider the thermal evolution of convecting bodies subject to 
strong viscoelastic tidal heating. We compare results first from simple models using a fixed Quality 
factor and Love number, and then for three different viscoelastic rheologies: the Maxwell body, the 
Standard Anelastic Solid, and the Burgers body. The SAS and Burgers models are shown to alter 
the potential for extreme tidal heating by introducing the possibility of new equilibria and multiple 
response peaks. We find that tidal heating tends to exceed radionuclide heating at periods below 
10-30 days, and exceed insolation only below 1-2 days. Extreme cases produce enough tidal heat to 
initiate global-scale partial melting, and an analysis of tidal limiting mechanisms such as advective 
cooling for carthlikc planets is discussed. To explore long term behaviors, we map equilibria points 
between convective heat loss and tidal heat input as functions of eccentricity. For the periods and 
magnitudes discussed, we show that tidal heating, if significant, is generally detrimental to the width 
of habitable zones. 

Subject headings: celestial mechanics — planetary systems — planets and satellites: general 



1. INTRODUCTION 

The discovery and study of planetary systems around 
different stars has revealed a rich diversity of orbital ar- 
chitectures, many of them not anticipated. Among the 
surprises is wide evidence for planet migration, orbital 
resonances, and the ubiquity of high orbital eccentrici- 
ties. These interactions may often lead to orbits that 
allow for durable nonzero eccentricities close to the star, 
especially for terrestrial mass planets. Our goal is to in- 
vestigate the range of tidal magnitudes that result from 
such orbital conditions. This paper examines the global 
temperature behavior of a simplified terrestrial planet 
during long-term extreme tidal heating, perhaps driven 
by mean motion or secular orbital resonances, using sev- 
eral different models of viscoelastic material response. 

We first present results across a range of orbit periods 
using a blackbody model with fixed material parameters. 
This method suggests that at eccentricities of around 0.1, 
tidal heating becomes globally significant generally at he- 
liocentric periods below 10-30 days, and large enough 
to drive global-scale partial melting at periods near 1- 
2 days. We map this range of extreme tidal activity in 
more detail using viscoelastic methods, and by modeling 
both temperature-dependent viscosity and melting. We 
investigate tidal work functions for three anelastic The- 



ological rock models: the Maxw ell, Standard Anelastic 
Solid (SAS), and Burgers bodies (fNowick fc Ber7vl ll972l: 
iSabadini et aLlll987t |Coorjer||2OO0). Each of these alter- 
native rock models have the potential to exhibit com- 
plex behaviors because their frictional work function is 
non-monotonic in temperature. Lastly, behaviors and 
extreme tidal equilibrium states are explored by linking 
heat input to a parameterized convection model. 

In this paper we use the term supertidal to re- 
fer to planets where tidal heating contributes signifi- 
cantly to the global heat budget. While the conditions 
needed for enduring extreme tides are expected to be 
rare, we are motivated by the rapid expansion of exo- 
lanet discoveries. Transiting s earches such as K epler 
Borucki et al1l200l. COR O T (IBorde et all [2001 and 
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MEarth ( Charbonneau et alj I2008T ) all have the poten- 
tial to discover a terrestrial planet of this supertidal 
type. A Hot Earth trapped in 2:1 resonance with a Hot 
Jupiter may also be detectable via careful transit timing 
(|Holman fc Murravi[200l : iMiller-Ricci et al.ll2008f) . how- 
ever no such timing patterns have yet been found. 

We focus here on hot planetary cases, where tidal 
heat competes with insolation fo r significance . This 
is as opposed to cold exomoons JScharf 200|), where 
it is easier for tidal heating to play a more dominant 
role, but detectability may be further off. Tidal heat 
may also p lay a role in enlarging the radii of some 
gas giants (jBodenheimer et all 120031 ). however we fo- 
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cus on rocky objects. Study of tidal cxoplanets is im- 
portant partly because eccentricities in the overall exo- 
planet population are higher t han expected (iButler et aTl 
120061 : iTremaine fe Zakamskall2004t iNamounil 12005( 1. e.g. 
GJ436b at e=0.15 ([Deming et alJ 12007( 1. suggesting the 
importance and prevalence of multiple planet interac- 
tions. 

Heating predictions for Hot Earths can become quite 
large, on the order of millions of terawatts. However 
global-scale partial melting of a mantle can begin with 
only a few tens of terawatts of added tidal heat. As on Io, 
widcscale melting lowers average viscosities and can lead 
to tidal shutdown. We discuss this shutdown behavior 
and how it can be altered by the high pressures in an 
Earth-mass body. 

Our methods for a Hot Earth follow thos e used to an- 
alyze the anomalously high heating of Io ([Smith et al.l 
[19791 ICarr et all [ 19981 IMcEwen et all 12000ft and to a 
lesser extent Europa dCassen et a!T 1979t iSauvres et al.l 
119831: lOiakangas fe Stevensonl Il989ft ~ and Enceladus 
( Porco et al.l 120061 1. These bodies have been the subject 
of intense a nalysis in terms of tidal heating via constant 
parameters (IRevnolds et al.1 Il980t lYoder fe Peald 119811 ; 
I Tackle v et all 2001ft. viscoelastic parameters and l ayere d 



def or mat io n (iSegatz et al.l 119881; iTakeuchi et al 
advection ([Moord 120011; iMonnereau fc Dubuffet 



1962 1. 



2002). 



and thermal-orbital coupling through orbital eccentricity 
(iHus smarm & Spohn 2004; Showman & Malhotra 1996; 
IPeale fc Ledl2002t IShowman fc Hanll2003ft . Here~we con- 
sider melting and tides as a function of temperature and 
forcing frequency, but do not explicitly model orbital 
feedback or include variations by layer, latitude or longi- 
tude. 

In section 2 wc discuss orbital configurations that can 
give rise to extreme tidal heating. In section 3 we present 
the model and results for a fixed Quality factor (Q) and 
Love number {k^) approach. In section 4 we discuss the 
derivation of complex rigidities, Love numbers, and work 
functions for three viscoelastic models and show results 
in section 5. In section 6 we discuss long term behaviors, 
tidal equilibria, and present bifurcation diagrams for the 
movement of tidal-convective equilibria points through 
changes in tidal forcing. Lastly, in section 7, we address 
issues such as advective heat transport, magma ocean 
production, and habitability. 

2. RESONANCE AND STABILITY 

The first condition for extreme tidal heating is proxim- 
ity to a massive host, providing a large change in gravity 
gradient between pericenter and apoccntcr. While moons 
often meet this criterion, only now have a large number 
of planets been detected in regions near stars where tidal 
heating becomes of geological concern. For extreme tides 
around a typical main sequence star, planets must be well 
inside the 88 day orbit of Mercury, but precisely within 
the 1-20 day range of Hot and Warm Jupiters. An ap- 
parent number density peak of Hot Jupiters exists near 
^4 day periods, and thus a mirrored peak in 2:1 reso- 
nant Hot Earths could occur near ^2 day orbits (unless 
proximity to the star causes better stability only further 
away). 

The second condition for tidal heating is an elliptical 
orbit maintained by some durable perturbation, such as 
the regular effect of Europa on Io. Orbital resonances, 



and in particular the 2:1 mean motion resonance, provide 
an effective mechanism to support long-term tidal heat 
production. Short-term tides from obliquity and non- 
synchronous spin damp quickly, leaving the body briefly 
warmed, but thereafter tidally quiescent. Nonzero ec- 
centricities can be more enduring, and can arise from 
mean motion resonances, secular perturbations, secular 
resonances, or recent interactions. 

2.1. Mean Motion Resonances 

Several mean motion resonances (M MR) have been ob- 
served in extrasolar planet syste ms (|Lee fc Peald 120021 : 
IGu et all 12001 Kiev et all [2005ft . For several reasons, 
the 2:1 MMR and to a lesser degree the 3:2 MMR are 
most favorable for exciting a terrestrial planet to ex- 
treme long-term tides. In a resonance {p + q) : p, the 
order q, determines how often orbital conjunctions are 
co- located (e.g. at pericenter). Frequent coherent con- 
junctions in the 2:1 and 3:2 cases increase resonance per- 
turbation strength relative to higher orders such as the 
3:1 or 5:2 resonances. The 2:1 MMR can be particularly 
stable because of its large physical gap in semi-major 
axis from other resonances, whereas the 4:5 and higher 
MMRs can become packed together in real space, leading 
to chaos and ejections. 

The large resonance widths of the 2:1 and 3:2 MMRs 
also favor longer durations for extreme tides. Since tidal 
heating will drive a system away from a state of perfect 
resonance, wider cases can induce greater heating prior 
to resonance breaking. 

The order of capture for a migrating perturber favors 
the presence of bodies in the 2:1 and 3:2 resonances. Cap- 
ture here is likely, and bodies must have either missed 
capture or have been otherwise scattered to reach and 
populate the resonances closer to the perturber. 

High perturber eccentricity can be both favorable and 
detrimental for extreme tidal forcing. High perturber 
eccentricity in the 2:1 resonance may increase the am- 
plitude of inner body eccentricity forcing. However, if 
pumping is too great, orbit crossings can lead to ejec- 
ti ons and clearing. 

iHaehighipour fc Raymond! ([2007( 1 and 

iRivera fc H aghighi pourl (|2007ft analyze the stability 
of terrestrial planets in or near resonances with short- 
period gas giants and known exoplanets and identify 
multiple orbits both inboard and outboard of the 
perturbing bodies as stable or quasi-stable. 

Once strong tides are established, they tend 
to force the system away from perfect resonance. 
iTerquem fc Papaloizoul ([2007ft analyze the dynamics of 
short period Hot Earths and find tides eventually break 
exact 2:1 resonances, but remain in near-resonance with 
nonzero eccentricities during later stage t idal dissipation. 
Disk turbulence may also break MMRs ([Lecoanet et al.l 
12009ft but leave bodies in a near-commensurable state 
vulnerable to ejections. 

2.2. Secular Perturbations 

Secular perturbations occur in two-planet systems and 
lead to an equilibrium forced eccentricity. In this case 
there is no resonance to be broken, so such forced ec- 
centri cities can proy i de a very long t erm source for 
tides (|Mardlingj|260l iZhou fc Lin|[200l) . The effective 
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strength of this perturbation against damping depends 
on the secular forcing timcscale (a function of the planets' 
sizes and locations) and the tidal dissipation timescale (a 
function of heating magnitude). 

Following iMardling fe Linl (|2004l ). the secular equilib- 
rium eccentricity for an interior Hot Earth of semi-major 
axis a sec perturbed by an outer Jupiter (not in 2:1 res- 
onance), is ~6% of the eccentricity of the perturber for 
apert/asec=10, and 19% for a per t/a sec ~30. Relativistic 
corrections can reduce these values somewhat further. 
Secular timescales are often of order 10,000 years, while 
moderate tidal timescales are often on the order of a few 
million years, suggesting that conditions can exist where 
modest geologically significant tidal heating is supported 
by secular perturbations alone. While extreme tides may 
damp forced eccentricities, temporarily cessations per- 
haps due to mantle melting will allow windows for secular 
forced eccentricities to be restored. 

Secular resonance occurs in multi-body systems. Un- 
stable secular resonances occur at fixed periods, which if 
a small planet crosses can lead to high eccentricities and 
ejection. In our solar system these unstable points lie 
due to Jupiter and Saturn near 0.5, 2, 12.5, and 17.5 AU 
(jMurrav fc Dermottl [20051 ) . Outside these singularities, 
secular resonance drives eccentricities of order 0.0-0.06. 
Unlike with MMR cases, there is no locking mechanism, 
so a planet under the influence of a secular resonance 
can migrate inwards due to tidal heating. If in one of the 
singularities, migration will be more rapid, leading to a 
brief episode of increased heat output. This effect will 
be magnified the closer to a star it occurs. Such secular 
singularity crossings are likely explanations for sudden 
bursts of eccentricity, but not of sustained extreme heat- 
ing. 

2.3. Capture 

Mean motion resonances in exoplanet systems arc 
thought to arise from various conv ergent mi grations in 
the late stages of planet formation (|Led 1200 4). Differen- 
tial inwards migration is com mon in early solar systems 
(iThommes fc Lissauerl I2005D and pla netary moon sys- 
tems (|Kleyill999HCanup fc Wardll2002f ) . and is a possible 
origin of the 1:2: 4 resonant Laplac e relation among the 
Jovian satellites (jPeale fc Ledl2002ft . Resonance capture 
occurs in convergent cases when eccentricities are below 
a critical threshold. Capture does not occur in divergent 
cases. For large initial eccentricities, capture becomes 
probabilistic, with a kn own, albeit complex , expre ssion 
for P C apture- Following iMurrav fc Dermottl (|2005f) . the 
critical eccentricity for internal 2:1 resonance capture 
in an Earth-Jupiter-Sol system is 0.152, and for a 3:2 
capture is 0.122. Threshold values can be higher for M 
dwarf and super-Jupiter cases. A terrestrial planet's ef- 
fective starting eccentricity is a function of both is secular 
forcing environment and its impact history. While tides 
are unlikely to cause much damping at 1AU during the 
chaotic growth phase, many early inner planets may have 
e < 0.15 and be candidates for capture. 

One scenario to form a supertidal Hot Earth is for 
a migrating future Hot Jupiter to capture and sweep 
an inner rocky planet along in its 2:1 inner resonance 
(|Yu fc Tremaindl200lD . For young solar systems, nebular 
disk torques are expected t o dominate and induce Typ e 
I and Type II migrations (|Thommes fc Lissauerl I2005D . 



After a gas giant clears a gap in its local region, Type 
II migration takes hold, and is considered the source of 
Hot Jupiters. Terrestrial planets may e xperience mod- 
est ra ndom-walk style Type I migration (jLaughlin et alJ 
12004ft . but remain generally in place, setting up a con- 
vergent pattern. Sweeping capture may collect multiple 
bodies from an inner solar system into the same reso- 
nance. The sweeping merger of multiple large embryos 
would lead towards a super-Earth mass with an initial 
hot start. 

It is favorable if the MMR reaches a rocky planet prior 
to a secular resonance that could pump up high eccen- 
tricities. Secular resonance positions depend on the over- 
all solar system configur ation and shift when p recession 
rates or masses change ([Nagasawa et al.1 120051 ). While 
many scenario geometries can prevent a secular singular- 
ity from disturbing a candidate 2:1 Hot Earth, one such 
case is a system with only one gas giant and thus a lack 
of unstable secular resonances. Si nce only 5% of sunlik e 
stars appear to have gas giants (jUdrv fc Santos! 12007ft . 
s uch cases may be comm on. 

lYu fc Tremahld" (|2001f ) and Eel (|2004ft analyze the res- 
onant sweeping of exoplanets in detail. As a giant outer 
planet continues to migrate, it progressively increases the 
trapped inner body's eccentricity, and later inclination. 
After migration by a factor of 4 in semi-major axis (for 
a 2:1 case), the inner body eccentricity is large enough 
to cause close encounters with the host star. Loss to the 
star can occur, or perhaps tidal damping and disk inter- 
actions can hold eccentricity to an equilibrium value as 
migration proceeds further. If the trapped Hot Earth is 
able to survive in resonance all the way down to short 
periods, it may do so with a large initial reservoir of ec- 
centricity to feed further tides. Undamped resonant mi- 
gration by a factor of 9 leads inner bodies to resonance 
release on nearly circular retrograde orbits. For these 
reasons 2:1 trapped Hot Earths may be more favorable 
at dimmer stars, where snow lines arc closer and requisite 
migration distances shorter. 

Alternately, orbital resonances can be traps where 
scattered planetcsimals congregate witho ut invoking the 
sweep ing capture mechanism above. iMandell et alJ 
(|2007ft demonstrate in numerical simulations of gas-disk 
induced migrations the formation of a variety of Hot 
Earths via scattering, often near the 2:1 resonance points 
of migrating Hot Jupiters, and with inner solar systems 
often cleared of further material. Denser inner gas disks 
appear correlated with having Hot Earths at the end of 
their 200 Myr simulations. These simulations did not 
include tidal damping or attempt to address long term 
stability. 

Outward migration due to the influence of a tidal bulge 
on a fast rotating primary is considered nec essary for the 
Galil ean system to remain in equilibrium (jPeale fc Led 
[2002] ). Analo gous behavior may occur for hot exoplan- 
ets. For the outward migration of a single body not in 
a resonance, eccentricity will simply increase. High e in 
moon systems is considered a signature of outward tidal 
evolution in older systems. When the joint migration 
of a 2:1 resonant system is dominated by dissipation in 
the primary, eccentricities will be limited by Q of the 
primary, which sets an effective limit on the d issipation 
in both secondaries (|de Pater fc Lissauerll2005ft . For our 
analysis we make no assumptions about tides in the star, 
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nor the star's rotation rate relative to its inner planets. 
We generally consider e as a free parameter and compute 
resulting tidal heat rates. 

Full treatment of tides requir es thermal-orbital cou- 
pling (Hussmann & Spohn 2004), which this paper does 
not model. In resonant systems this coupling is complex, 
and extreme friction can lead to evolution away from ex- 
act resonance and possible resonance breaking. In non- 
resonant systems tides simply lead to circularization and 
semi-major axis drift. 

2.4. Circularization 

While ongoing perturbations are favorable to super- 
tidal conditions, they are not necessary. Circularization 
timescales may still be of t he order 0.1-10 Gyr for short 
period Earth-mass planets. iJackson et alj (|2008[ ) also ex- 
amine the tidal heating of non-resonant terrestrial exo- 
planets, and discuss how tidal orbital migration further 
lengthens circularization times. 

Often circularization times are reported without regard 
to a planet's internal response, leading to results that 
may be too short. For example a Hot Earth in a 2 day 
orbit of a 1 solar mass host with 1 x 10 4 terawatts (TW) 
of dissipation at e = 0.1 will have a 5.4 billion year damp- 
ing timescale, while at lxlO 6 TW circularization takes 
only 54 million years. The higher dissipation rate how- 
ever may be unsustainable for a rocky body due to global 
onset melting (or more generally viscoelastic de-tuning) . 
Instead, a migrating planet may briefly experience a few 
million TW in tides, then rapidly warm to a small partial 
melt fraction. Decreasing viscosities thereafter produce 
only a few thousand TW of tidal heat, which may then 
persist for a solar system lifetime even without support 
from a resonance. This phenomenon is discussed further 
in section O It is important when calculating circular- 
ization timescales using a fixed Quality factor Q to take 
into consideration whether the planet itself can geologi- 
cally produce the subsequent rate of dissipation. 

Countless complications to the arguments in this sec- 
tion exist: High exoplanet eccentricities could in part 
be caused by unseen high inclin ation perturbers, as oc - 
curs via the Kozai mechanism (Takcda fc Rasiol I2005D . 
Nonzero secular precessions and inclinations mix with 
mean motions to obtain a variety of higher order reso- 
nances. At short periods, a star's oblatcncss, tidal bulge, 
atmosphere, and relativistic effects may also influence 
stability. 

Observations will ultimately decide the matter. Over- 
all we consider it likely enough that some terrestrial plan- 
ets can be swept or scattered into resonances by migrat- 
ing Hot Jupitcrs, or may otherwise have their eccentric- 
ities sustained at nonzero values for geologically signifi- 
cant times, to move forward and consider the tidal heat 
magnitudes that then result. 

3. FIXED Q TIDAL MODEL 

Tidal heating is modeled in many ways, but 
the starting point is the global heat generation 
rate E ud ai (iPeale fe Cassenl 119781: IPeale et all 119791: 
IShowman fe Malhotral 119961 ). For a homogeneous spin- 
synchronous body whose stiffness and viscous dissipation 
are both assumed to be constant and uniform, the global 
tidal heat ra te can be expressed follow ing the detailed 
derivation in lMurrav fc Dermottl (|2005l ). 



E, 



21 k 2 GM^„R 5 sec ne 2 
2 Q a e 



(1) 



Eudai'- Tidal heat production rate, watts 
a: Semi-major axis 
e: Eccentricity 

G: The gravitational constant 

Mp r i : Mass of the primary 2 

R S ec- Radius of the secondary 

k 2 : Second-order Love number of the secondary 

Q: Quality factor of the secondary 

The masses, radii, and orbital parameters a and e for 
exoplanets have the potential to be directly measured. 
In contrast, Q and k 2 must be inferred, and represent 
in only two (albeit convenient) lumped scalar terms all 
the complex internal material properties of a deform- 
ing planet. After the high power terms in equation [1] 
are known, almost all uncertainty in the tidal heating 
comes from k 2 and Q. As discussed in section [4j these 
terms together are actually a function of both forcing fre- 
quency and temperature. But first we can gain insight 
by demonstrating results of equation [1] simply by picking 
commonly invoked values for k 2 and Q. 

The Quality Factor Q is an inverse damping term char- 
acterizing energy lost to friction, and can be variously 
defined for a ma terial, mode, process, or c ompo und pro- 
cess. Following lO'Connell fe Budianskvl l|1978f) . the Q 
for an oscillator may be written: 



Q{u) = 2tt 



(2E ave (u)) 
AE(lu) 



(2) 



Where E ave is the average stored energy, u> the fre- 
quency, and AE the energy lost per cycle. This differs 
slightly from definitions of Q using E max , the max imum 
stored energy per cycle, but following Bland (1960) equa- 
tion's best compatible with the formal viscoelastic def- 
inition of Q we discuss in section [H equation [T71 Note 
that Q will differ depending on the frequency for which 
it is measured, as well as on the temperature and micro- 
dynamic properties of a material. 

F or the whole Earth, a Quality factor of 12 to 
34 (jYoderi [l995l iMurrav fc Dermottl I2005T ) has been 
measured for the overall tidal p rocess based on th e 
expansion of the moon's orbit (|Dickev et alj I1994T) . 
These values however include all damping effects of 
the lunar tide on the Earth, including the sloshing 
of Ear th's oceans, and v aried continental configura- 
tions. IGoldreich fc Soteri (|1966D estimate Q for Mer- 
cury, Venus, and the Moo n at Q m erc<\ 9 0i Qv <17, and 
10<Q r „<150 respectively. lLainev et al.l (|2007f ) recently 
estimate Q mars =79.91 ( ±0.69) from the mot i on of Pho- 
bos. Data summarized in lKarato fc S petzlcr (199fJ) sug- 
gest Earth's lower mantle Q at tidal periods from 1-10 
days is in the range 50-200. Widescale partial melting 
may also drive ocean-free 1Mb Q's closer to 10. For an 
evolved dry Earth-mass planet, we adopt a solid body 

2 Because subscripts such as M s and M p can be confused be- 
tween planet, satellite, sun, primary and secondary, we adopt three 
letter subscripts and use primary to refer to the central mass, 
whether it is a star or planet, and consider tides generated inside 
an orbiting secondary. 
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baseline of Q=50. Extending such results to Q=10 or 
Q=100 is straightforward given Q is linear in cquation[TJ 

The Love number fc, is in essence a compliance term, 
characterizing the deformation response of a planet to 
stress. As tides excite the 1=2 prolate harmonic, we are 
concerned only with the &2 term. A of represents 
a perfectly rigid bo dy, and 3/2 a perfect fluid. Earth 
has a value of 0.299 (|Yoderill995l ). Through self-gravity, 
fc2 depends strongly on an object's size, with Mars at 
0.14, and Jupiter still far fr om the fluid limit at ^2=0. 38 
(|Gavrilov fc Zharkovlll977f) . Small strength dominated 
moons such as Phobos and Amalthca have Love numbers 
near zero. Thus for most terrestrial planets there is less 
than a factor of 2 uncertainty in fc'2 . 

As we proceed it is important to keep in mind errors in- 
curred by our model choice. Equation 1 is based on linear 
tidal theory, and we are applying it to cases where dissi- 
pation is not small. Exoplanet stains in the region of in- 
terest for tidal heating range from 10~ 8 to as high as 10 -4 
or 10~ 3 in extreme cases, depending on eccentricity, while 
nonlinearities may become important at strains greater 
than 10 -6 . Corrections involving higher order terms and 
small angle approximations such as those discussed by 
lEfroimskv fc Williams! jUS) for tides in a primary, are 
more likely to become non-negligible for extreme tidal 
exoplanets. These corrections however may themselves 
be small compared to the greater error of not including 
planet-wide partial melting. 

Values of &2=0.3 and Q=50 at lMg allow us a simple 
baseline to highlight other features of the tidal heating 
parameter space. 

3.1. Energy Balance 

To judge the geologic importance of tidal heating, 
we compare against other heat sources, specifically in- 
ternal radiogenic heating and insolation from the host 
star. We estimate surface temperatures based on a black- 
body assumption and a summation of primary energy 
sources: insolation, radionuclides, accretion (or gravita- 
tional) heat, and tides: 

Efotal — -^insol ^radio ~t~ Eg rav ~f~ -^tidal (3) 

The blackbody surface temperature of such a model 
planet is found from the Stefan-Boltzmann Law 

Etotal = Euackbody = 47T R 2 sec O - s be r T^ ur j (4) 

where a s b is the Stefan-Boltzmann constant, and emis- 
sivity e r is close to unity (0.9 in the infrared). By con- 
sidering planets as blackbodies we implicitly assume thin 
atmospheres and no significant greenhouse effects. 

We have tested models both with and without diurnal 
variations, since supertidal planets are likely to be spin 
locked to their host stars. Hcliosynchronous spin lock- 
ing may establish a strong boundary condition contrast 
for the mantle, of 1000K or more between dayside and 
nightsidc, if held stationary long enough. Such contrasts 
may alter convection plume patterns and induce true po- 
lar wander of vigorous elevated limb plumes back to the 
equator. 

Equating E to tai to the outgoing radiation is an assump- 
tion of steady state. The time it takes a given supertidal 



planet to reach equilibrium is discussed in section but 
is generally on the order of a few million years. 

We scale insolation from the host star luminosity, L star 
in watts. Correcting for albedo A, the total insolation 
heat rate of a planet is: 

£ mso! = ML)(l-i)^j (5) 

This radiation its rapidly re-radiated at thermal wave- 
lengths, but in so doing it adds to the outward black- 
body flux, establishing a surface boundary condition for 
the planet's internal gcothcrm. 

Estimates fo r Earth's internal h eat flux range from 
30 to 46 TW (|Jaupart et al.l I2007T) . Measurements are 
made by sampling the geotherm in deep sea cores and 
equilibrated drill holes, as well as satellite IR photome- 
try of the entire Earth. Much of the uncertainty comes 
from estimating the flux from volcanic point sources and 
through ocean crust hydrothermal circulation. About 
half of this heat is thought to come from radionuclides. 
The rest comes from several sources, including leftover 
accretion heating (gravitational potential heat trapped 
by burial), chemical phase change heat, primordial grav- 
itational heat of the core differentiation event, continuing 
latent and gravitational heat release from the inner core 
crystallization, heat from the moon forming impact, and 
the pulse of heat provided by extinct short lived nuclides 
such as 26 Al. In the deep Earth, tidal heating is negligi- 
ble. 

Although we are mostly concerned with ~1Me exo- 
planets, for super-Earths we scale radiogenic and residual 
accretion heat by the total planet mass relative to 1Mb 
(assuming a high silicate mass fraction). Using chon- 
dritic concentrations of 235 U, 238 U, 232 Th, and 40 K, we 
find radionuclide heating just under 9 times greater for 
an early Earth. Ac tual mantle nuclide concen trations are 
highly uncertain (jMcDonough fc Sunl 119951 ). However 
for tides to dominate an evolved planet they may need 
only exceed ~40TW, but to dominate a younger planet 
just at the end of its migration phase, they may need 
to exceed ~400TW. For radionuclide poor Ice Earths, 
smaller amounts of tidal heat will have greater conse- 
quence. We do not explore any range of core mass frac- 
tions or the possibility of radionuclide concentrations sig- 
nificantly different from our solar nebula. This in effect 
assumes a turbulent interstellar medium well mixed by 
random su pernovae, and is general ly supported by ob- 
servations (jElmegreen fc Scaloll2004f ) . However, the time 
since the nearest supernovae can vary the concentration 
of 26 Al and hence the initial pulse of heat to a planet. 
The occurrence or absence of an early giant impact might 
have a similar effect. 

3.2. Fixed Q Results 

Figure Q] compares the ratios of tidal heat to insolation 
and radiogenic heat for hypothetical Hot Earths (desig- 
nated by the suffix x) trapped in 2:1 resonances with 
known short period exoplanets as taken from the ex- 
oplanet. eu database of Jean Schneider. Scatter of the 
points is due to the varied luminosity of certain stars, 
with higher outliers being M dwarf hosts. A sufficient 
atmosphere is assumed to transport heat evenly to the 
nightside. 
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Fig. 1. — Comparison of tidal heat to radiogenic and insolation 
heat for hypothetical Hot Earths in 2:1 resonance with known ex- 
oplanets. Body parameters including the modern radiogenic heat 
rate (~19TW) for Earth are used. Even given a conservatively 
strong forcing of e=0.1, tides only reach parity with insolation 
for extreme objects inside 2 day periods, where global-scale par- 
tial melting may invalidate the fixed Q assumption. lAfg, Q=50, 
fc 2 =0.3, A=0.3. 

Tidal heat dominates over radiogenic heat at periods 
under about 30 days, while it only dominates insolation 
below ~2 days. This defines our range of interest for 
tides. When tides dominate radionuclides, they dom- 
inate the internal behavior of the planet. For super- 
Earths, large hosts, and uncommonly high eccentricities, 
the range of internal relevance shifts outwards, to typi- 
cally no more than ~80 day periods. When tides dom- 
inate over insolation, they can begin to be expressed in 
the observable temperature of the planet. 

These data show it is relatively easy for tides in short 
period objects to have a strong influence on internal ac- 
tivities, but it is very difficult for tides to matter for 
overall surface temperature. This is a consequence of 
the fact that to get strong tides in the first place, you 
have to be very close to a star, where Emsoi ~lxl0 6 
to 1 x 10 s TW. At dimmer stars, tides can have greater 
surface significance. 

Figure shows tidal relevance with respect to surface 
temperature. Again, the point is that as an observable, 
the surface temperature of even extreme tidal objects is 
generally overwhelmingly dominated by insolation. Sur- 
face temperatures with and without tides are nearly 
equivalent except for the very nearest bodies around 1 
day orbits, where the tidal signal reliably separates from 
the insolation signal. The plotted error bars represent an 
uncertainty range in albedo of (dark ash filled skies) to 
0.8 (bright solid clouds), with 0.3 (Earthlike) being the 
centermark. For anything but the most extreme cases or 
dimmest stars, uncertainty in albedo will prevent over- 
all surface temperature from being used as a means to 
observationally identify planets as supertidal. Beyond 2 
day orbits, the tidal contribution to blackbody temper- 
ature is negligible. Within 1 day periods, we argue that 
the tidal signal is erroneously high, due to the limitations 
of the fixed Q model with no melting or viscoelasticity. 

Even extensive volcanism induced by tides may have 
little impact for close planets of G class hosts, since the 
stars alone can maintain high enough temperatures to 
liquefy solid material. For brighter A and F stars tides 
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Fig. 2. — Blackbody surface temperatures for hypothetical su- 
pertidal planets below 4 day orbital periods. Simulated data using 
the known catalog of exoplanets as 2:1 perturbers, with scaled stel- 
lar host masses and luminosities. Black marks with error bars: No 
Tides (e=0) with albedo uncertainty from A=[0.0, 0.3(nominal), 
0.8]. M dwarf hosts form the low outliers. Squares: Fixed Pa- 
rameter Tides results using 1Mb, 6=0.1, Q=50, A)2=0.3. Crosses: 
Viscoelastic Maxwell rheology (more conservative estimates). Note 
tidal surface enhancements are generally lost in albedo uncertainty 
(even at this extreme e) except at dim stars and periods below 1 
day. 



should matter even less. The zero pressure melting point 
of many silicates lies between 1400 to 2100 K. There- 
fore from Figure [2 insolation driven magma oceans are 
plausible, pending atmospheric effects, inside of 3 day G 
class orbits. The depth of such oceans will depend on 
the subsurface heat flux, ocean convection rate, compo- 
sition, and pressure dependence of the solidus. For large 
planets such as Earth, pressure rises rapidly with depth 
and prevents melting even at very high temperatures. 
Tidal heat, though minor compared to insolation, may 
contribute to increasing otherwise shallow magma ocean 
depths. The behavior and tidal sloshing of these insola- 
tion established magma oceans, and how they contribute 
to the global Q value is a subject of future work. 

For M dwarf stars and nonluminous hosts, tides mat- 
ter more. The reduction in tides due to a smaller host 
mass is easily outweighed by reduced insolation. It may 
be possible to observationally detect a few degrees tidal 
temperature enhancement at these dimmest stars. Note 
the data points labeled GJ581bx, GJ436bx and GJ876dx 
in Figure [2] A broader axiom is that tides matter most 
in colder environments. A relatively low 20 TW tidal 
output at Io is strongly expressed due to the absence of 
strong radionuclide heat and comparatively weak inso- 
lation. Although this paper focuses on planets, we do 
expect tidal heating is of greatest consequence, creating 
new habitable regions, in cold exomoon systems analo- 
gous to Io, Europa, and Enceladus. 

For spin-locked objects, even if only a small fraction 
(~0.1-1%) of insolation is transported to the nightside, 
large tidal contributions still appear negligible in most 
cases. If zero insolation is transported to the nightside, 
then opportunities do exist to detect the small tidal tem- 
perature elevations of 1-5 degrees that occur in orbits 
from 2-5 days, however this signal is still weak and such 
objects may be rare. 

Figure [2] suggests that to observationally confirm ex- 
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treme tidal heating near stars one should look for spec- 
troscopic evidence, such as magma hotspots or a chem- 
ical signature of widespread volcanism. Widespread dc- 
volatilization of a planet may produce extended C0 2 
and sulfur compound absorption lines in observed spec- 
tra, L. Kaltcneggcr et al. (2010, in preparation). Light 
from hotspots and lava lakes may be polarized as on Io 
(jVeeder et all 119941) and shifted to the near IR, aiding 
in detection. Clouds and haze may obscure hotspots, 
magma lakes, or even magma oceans, but if cloud breaks 
exist, then such observations are possible. In the mean- 
time, eccentricity and resonances are the best indicators 
of tidal heating. Chemical and hotspot follow up obser- 
vations can help in bounding global Q values. 

The most Io-like exoworlds will be those where in- 
solation is too weak to melt the surface on its own, 
but tides are great enough to drive volcanism in excess 
of background rates. Consider that E t idai / Ei nso i and 

Etidal/Kadio are -0.08, -77 for Io, and -0.01, -10 for 
Europa. Thus dominance over insolation is unnecessary 
for Io-like supertidal properties. While moderate tides 
may not noticeably change a planet's average blackbody 
surface temperature, they can slow secular cooling and 
radically alter layer structure, plate tectonics, and the 
state of volcanism. Moderate tides, of a few hundred 
TW, are also less likely to break resonances or result in 
rapid circularization. 

The main point from this exercise that motivates our 
further work is the magnitude of tidal heating predicted 
by the fixed Q and k 2 method. In Figure Q] the tidal 
heating of the shortest period hypothetical planets rises 
to over a million times the radiogenic rate: tens of mil- 
lions of terawatts. The question we seek to explore is 
whether or not such extreme answers are possible. Can 
an Earthlike planet experience millions of terawatts in 
tidal output for any length of time? A number of physical 
effects conspire against this extreme answer: viscoelastic- 
ity, global-scale partial melting, convective feedback, and 
thermal-orbital coupling. Next in this paper we address 
viscoelasticity and melting. 

4. VISCOELASTICITY 

Using equation Q] to calculate global tidal heat is use- 
ful for estimates, however it ignores the frequency de- 
pendence of a material's response to loading. Q and k 2 
are n either constant nor entire l y independent par ame- 
ters, (|Kaulalll964t IZschau|[l97l1Segatz et"aL|[l988l ). In 
its most general form for a viscoelastic body, the ratio 
k 2 /Q is replaced by the imaginary part of the complex 
Love number — Jmffe), which characterizes the mate- 
rial's viscous phase lag. 



Eudai = -Im{k 2 ) — — = 



(6) 



This formula still assumes a homogeneous body. A 
complete calculation of tides would consid er variations by 
layer s using a propagator matrix method (jTakeuchi et alJ 
119621 ) as well as the full three dimensional stress and 
strain tensors to compute tides as a function of latitud e 
and longitude (jPeale fe Cassenlll978l ; ISegatz et al.lll988l ). 
However equation [6] is effective in seeking estimates and 
extrema of a globally averaged behavior. 
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Fig. 3. — The four different anelastic models discussed in this 
paper, and their corresponding notation. 

In general, Im(k 2 ) will have a response peak at a ma- 
terial's characteristic natural frequency, allowing bodies 
to tidally resonate. To compute Im(k 2 ) for any rheology, 
we begin from a definition of the Love number: 



1 



21 + p, 

Where jl is known as a body's effective rigidity 



V = 



19/j, 



2pgR a 



(7) 



(8) 



and [i is a material rigidity. Effective rigidity is a non- 
dimensional ratio that compares elastic forces in the nu- 
merator to gravitational forces in the denominator. For 
ft >> 1 the body is strength dominated, while for fl << 1 
it behaves like a gravitating fluid. Often the substitution 



P = pgRs 



(9) 



is made to emphasize this point by forming a gravi- 
tational stiffness analogous to the material modulus /i. 
Since tides excite a planet in shear (orthogonal tension 
and compression), for we use the Shear modulus G 
(not the Bulk or Young's modulus K or E). Typical val- 
ues of G are 90 GPa for ( undamaged) rocky material , and 
4 GPa for icy material (jGoldsbv fc Kohlstedti[200ll ). In 
comparison, the effective gravitational rigidity is /J=250 
GPa for Earth and 8 GPa for Io. 

A material may behave as any combination of springs 
and dampers that model viscous creep and elastic re- 
bound. Mathematically, a spring damper model is rep- 
resented by forgoing a scalar modulus fi for a complex 
value M*, whose real and imaginary parts define the en- 
ergy storage and energy loss aspects of the system. These 
components Mi and M 2 are computed from a model's 
constitutive equation. Substituting the complex form of 
the stiffness 



fx(oj) = M*(w) 



Mi(w) + «M 2 (w) 



(10) 



into equation [8l then substituting this complex from 
of /2 into equation [7l we derive the complex form of k 2 
and extract the imaginary component Im(k 2 ) for use in 
equation [6] 

4.1. Anelastic Models 

The four basic models and parameters for the rhcolo- 
gies considered for this paper arc shown in Figure [3J The 
Maxwell body, commonly used because of its simplicity, 
considers mantle rock as a spring-dashpot series, with an 
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TABLE 1 

Viscoelastic Formulae 





Maxwell 


SAS 


Burgers* 


T 


V 

M 


ri&J 


T A = VA 5J, T B = $L 
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Mri 2 tu 2 


1 (v^fSJ 3 J u 


w (,W — VA" Ji ^2) 




Jr j r +j r (vuJ u &j) 2 


C 2 +u> 2 C 2 


M 2 


M 2 riu 


Vui&J 2 


^(C 2 +vaSJ^ 2 C 1 ) 


M' 2 +v 2 ^ 2 


J 2 +(V^Ju&J) 2 


C 2 +u> 2 C 2 


-7m(fc 2 ) 




576J 2 Pv^ 


57f3u(C 2 +SJu 2 C 1 v A ) 


4,p(l+[(l+(19M/20))*r)*u*/M*]) 


(19+2l3J R Y' ! +S.r- ! v 2 <<> 2 (19+20 Jm) 2 


361u 2 +76f3uj 2 C 1 +4(3 2 u, 2 C'{ + (2l3C2-196Juj 2 vA) 2 


Q 


V*jJ 
~M~ 


J u +SJ+J u (uvSJ) 2 


u(C 1 -vaWC 2 ) 


VtijSJ 2 


C 2 +vaSJ^ 2 C 1 




* d = 5J + SJHA + I 
VB M B 


* C 2 = J- - 2A«a, 2 
VB M B 



Note. — For the SAS model J u represents the initial elastic compliance, 57 is the additional later compliance due to creep, and 
Jr = J u + 5J is the relaxed compliance following creep. Burgers formulas are expressed in terms of the two lumped terms C\ (which is an 
effective compliance) and C 2 . 



instantaneous clastic response, followed by viscous yield- 
ing. It is ultimately fluid. 

The Maxwell model is useful but inco mplete 
( Oiakangas fc Stevenson! 119891 ) . As discuss ed by iZenerl 



( 19411) and observed in the laboratory (IPostl 1977 



Smith fc Carpenter! 119871 : iCooperl 120021 Uackson et ail 
2000f ) . real polycrystalline materials exhibit a wider 



range of relaxation mechanisms. 

A parallel spring-dashpot pair is known as the Voigt- 
Kelvin model. Here viscous relaxation is ultimately lim- 
ited by the spring. While the Voigt-Kelvin model is in- 
structive as a subcompone nt of other models and has 
been applied to lunar tides (|Van Arsdaldll98i| ). we ulti- 
mately find it poorly suited to short period cases. 

A three parameter model, known as either the Stan- 
dard Anelastic Solid, (SAS) or Standard Linear Solid, 
has features of both the Maxwell and Voigt-Kelvin prim- 
itives: instantaneous elastic response, followed by strain 
limited relaxation. It will not take a permanent set. All 
deformation is recovered when a load is removed. Either 
of the two ways to arrange two springs and one damper in 
a ser ies-parallel combinatio n are mathematically equiva- 
lent (jNowick fc Berrvj[l972l ). 

A four parameter model, or Burgers body, allows the 
modeling of transient molecular creep behavior in min- 
erals. It can exhibit transient creep, recovery, and take 
on a permanent set, modeling a broad range of materi- 
als. The Burgers or SAS models may both be reduced to 
the Maxwell or Voigt-Kelvin models through appropriate 
selection of parameters. 

The Burgers body is useful in modeling the phe- 
nomenon of grain-boundary slip. The Maxwell ele- 
ment within the Burgers body represents classical dif- 
fusion creep, where non-recoverable creep motion occurs 
through void migration inside of grains. Grain bound- 
ary slip occurs on a shorter relaxation timescale and 
is recoverable, as represented by the Voigt-Kelvin el- 
ement. Postglacial rebound studies in particular have 
suggested that the Burgers body is a m ore appropriate 
model of the Earth than a M axwell body ( Sabadi ni et alj 
Il987t iFaul fc Jackson! I2005h . As we apply them to exo- 
planct tidal models, both the SAS and Burgers mod- 
els reveal susceptibility modes not found in a Maxwell 
approach. We note that as a method of extrapolation, 
spring-dashpot models have the advantage of a clear re- 



lationship to underlying defect microdynamics, but the 
disadvantage of historically poor correlation with large 
complex inhomogeneous system s such as Earth's mantle . 

We follow the notation of iNowick fc Berrvi (|1972f ) 
where J's denote compliances in Pa^ 1 , and M's denote 
stiffnesses in Pa. Series springs arc mainly expressed by 
stiffnesses and parallel springs by compliances. For the 
SAS model J u is the instantaneous compliance to an ap- 
plied load, and 6J is the additional compliance during 
creep, known as either the creep defect or compliance 
defect. 

Using the geometry and parameters of each model from 
Figure 02 then solving for the overall stress a and strain 
e leads to the following constitutive relations: 



Ma + f]ij = Mr]i 
Ja = e + rjJe 
(J u + SJ)a + i]6J J u a = e + rfSJe 



(11) 
(12) 

(13) 



Maxwell 
Voigt-Kelvin 
SAS 
Burgers 

^ + (sj + 5 J^ + ^-)a+^a = e + VA SJe (14) 
rjB \ vb M B J M B 

While a tidal distortion of a planet is best expressed as 
an applied strain, solutions were also investigated for ap- 
plied stress problems. The main difference between these 
methods is that the Maxwell and Voigt-Kelvin models 
switch behaviors. In a Maxwell model there is nothing 
to prevent infinite viscous work under an applied stress, 
as no spring limits the dashpot's travel. Under an ap- 
plied strain, work is inherently finite. The reverse is true 
in a Voigt-Kelvin model, where rapid forced strain cycles 
can drive the dashpot to approach infinite work. Mainly 
for this reason we find the Voigt-Kelvin model unsuited 
to tidal cases, and the most insight is found from com- 
paring only the Maxwell, SAS and Burgers models under 
a cyclic applied strain. 

The generalized equations for cyclic applied strain forc- 
ing are: 



a(uj) = (di + icr 2 )e 1 ' 



(15) 
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e(w) = e e^ (16) 

Taking the first two time derivatives of equations 1151 
and!16[ then substituting these derivatives into the con- 
stitutive equations, leads to the complex stiffnesses Mi 
and Mi as they are reported in Table [1] 

Table [T] also reports the quality factor, relaxation 
timescalc r, and imffe) formulae of each model. The 
Quality factor Q{r\,Lo) is found as the ratio of the com- 
pliances (jO'Connell fc Budianskvlll978| ) , and is the same 
value under either applied stress or strain. 



Ji{y,Qj) = Mi (77, a;) 
J 2 (?7,w) M 2 {t] 1 uj) 



(17) 



Work per cycle per unit volume for a generic object 
can be found from the integral of stress times strain. In 
a sinusoidal uniaxial applied strain case: 



W(r)(T), u) =irM 2 ( V (T),u)el 



(18) 



This formula is nearly equivalent to equation [5] once 
tidal strain is computed from the orbit, however it omits 
key ingredients, specifically: self-gravity, a shape factor 
of 2/5 for a sphere, and a factor of 7/3 for the combina- 
tion of what are termed the radial and librational tides 
(arising from spin synchronization in an elliptical orbit). 
Self-gravity in effect acts like an additional spring in all 
of the diagrams in Figure [31 When a bulge is raised on a 
planet's surface, both gravity and material rigidity act to 
pull the bulge back down. To properly include self grav- 
ity, we use equation [6] and express complex compliances 
in terms of 7ro(fc 2 ) following equation [7] and the method 
outlined at the head of this section. These results for 
ImQzi) are given in Table [TJ 

Using equation [6l W(r)(T),u) = E t idai is entered into 
a global energy balance along with radionuclides, convec- 
tion, melting, and a viscosity model, then propagated in 
time to observe the temperature behavior of the planet. 



4.2. Viscosity Model 
For a viscosity model we use an Arrenhius relation 



(19) 



where R is the universal gas constant, E* an activation 
energy, and r\ a defining viscosity. 

From isostatic rebound studie s of Canada and Scan- 
dinavia (jMitrovica fc Fortell2004n . Earth's viscosity is of 
the order ~ 4 x 10 20 to ~ 1 x 10 23 Pa-s from the upper 
to lower mantle respectively. Results can be highly sen- 
sitive to this parameterization. We choose r\ to match 
a selected setpoint r/ se t, either 1 x 10 22 Pa-s at 1000K 
(a weaker or wetter case), or 1 x 10 24 Pa-s at 1000K (a 
stronger or more devolatilized case). Higher r] set choices 
in general attenuate the global response. Liquid phase 
viscosit y is modeled with magm as near 1000 Pa-s at 
2000K (jMcBirnev fc Murasdfl98l . 

The energy barrier for some combination of creating 
and moving lattice defects in creep is characterized by 
E* . Values relevant for dislocation creep at the stresses 
and pres sures of Earth's mantle are from 300 to 40 
kJmor 1 (lAshbv fc Verralll[l977l : Ifischer fc Spohnlll990fl . 
At very low stresses the rhcology may switch to diffusion 



creep, eit her Coble or Nabarro- Herring based on the tem- 
perature (jFrost fc Ashbvl[r982| ). however we assume the 
stress state for supertidal planets remains high enough 
(1 GPa or more) to neglect these creep mechanism tran- 
sitions. In addition we neglect pressure variation of E* 
and changes due to varied concentrations of dissolved 
volatiles. 

Baseline parameters are summarized in Table [2j Given 
77 (T) it is possible to replace the constant terms r\ in Table 
[jand obtain tidal work as W(r](T), ui). 

While not modeled in detail here, ice viscositie s range 
from 1 x 10 11 to 1 x 10 14 Pa-s (IPoirier et al.|[l98l . With 
a 4 GPa rigidity, this produces Maxwell times from a 
few minutes to a year. Thus Ice Earths are strong candi- 
dates for a peak viscoelastic tidal response. In addition 
to warming icy moons outside the snowline, tides may 
play a unique role in the degradation of ice mantles on 
large migrated inner solar system objects, as a way to 
emplace heat beneath thick high albedo cloud decks. 

Earth's bulk resonant Maxwell timescale r (see Table 
[1]) is a few thousand years, fairly unresponsive to a tidal 
period of a few days, and with an attenuating im(fc 2 ) 
factor of^9xl0 -6 . Io at 1.7 days is suggested to 
be nearer its Maxwell response peak, meaning a mantle 
viscosity ~ 1 x 10 16 Pa-s, perhaps due to temperature, 
partial melting, and composition. Earthlikc bodies may 
require a heat pulse trigger, or pre-softening, to reach 
such a tidally responsive Io-like state. Impacts, impact 
induced librations, and impact induced obliquity tides 
have been sugge sted as triggers for ti dal heating in Sat- 
urnian satellites (|Castillo et al.ll2006al lbl). Post -migration 
planets may also experience tidal heat while still hot from 
accretion and 26 Al. Once well coupled, tides have the po- 
tential to keep mantles warm thereafter. Alternatively, 
such triggers may not be needed to strongly couple tides 
in the SAS and Burgers models, where response peaks 
occur at lower temperatures due to the term SJ. 

4.3. Melting Model 

A description of silicate melting allows us to resolve 
both the rapid increase in convective vigor and the de- 
coupling of tides that simultaneously occur when vis- 
cosity and shear modulus decrease. P arametr i c mod - 
els of melting for Io are p resented by iMoord (|2003bl ) 
and iFischer fc Spohnl (1 19901) based on laboratory exper- 
iments by iBerckhemer et al.l (|1982[ ). These models vari- 
ously represent the essential feature of a breakdown tem- 
perature: at some point in the partial melting (or crys- 
tallization) process, a material switches from being best 
described as a solid matrix with fluid pores, to a fluid 
bath with isolated floating crystals grains. When grains 
loose contact with one another, the material looses shear 
strength and switches to the viscous properties of the 
fluid. 

We follow the IMoord (|2003bl ) model melting values of 
T sol = 1600K, Tug = 2000K, and T brkdwn = 1760K (40% 
melt), 1800K (50%melt), or 1840K (60% melt), along 
with a mantle specific heat C v — 1260 J/kg-K, and a 
heat of fusion Hf = 500,000 J/kg. Shear modulus G is 
constant up to the solidus then begins to drop according 
to an empirical Arrhenius law 



G = Goe 1 - 



-25) 



(20) 
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TABLE 2 

Baseline Material Parameters 



Series Shear Modulus 


M, M B = 


5 x 10 iu 


Pa 


Creep Compliance 


SJ = 


4 x 10~ 12 Pa" 1 


Defining Viscosity 


Vset = 


1 x 10 22 


Pa-s 


Burgers Parallel Viscosity 


VBset = 


2 x 10 20 


Pas 


Activation Energy 


E', Eg = 


300 kJmol" 1 


Solidus Temperature 


Tad — 


1600 K 




Liquidus Temperature 


Tu q = 


2000 K 




Breakdown Temperature 


Tbrkdwn 


1800 K 





Note. — Used for reported heat rates unless otherwise specified. 

with an activation temperature 40000K. Upon break- 
down, the shear modulus drops to zero, and the viscos- 
ity immediately becomes that of the liquid. Between the 
solidus and the breakdown temperature, viscosity drops 
off by a power law 

V = VsoUde- 40 * (21) 

where x is the percent melt fraction and r] S oUd is the 
viscosity at the solidus temperature. Moore also presents 
a stiffcr material where the empirical coefficient -40 is 
replaced by -10. 

The primary shortcoming of our melting approach is 
applying it to the planet as a whole. This model assumes 
a uniform adiabatic mantle temperature, a single solidus 
temperature, uniform global onset melting, and ignores 
compositional, pressure, eutectic, and hydration effects 
on the solidus. This model does allow us to capture the 
gross behavior of tidal shutdown at some small degree of 
partial melt, without telling us how the melt distribution 
affects the exact tidal result. In section [5] we discuss the 
role of local partial melting as well as the role of heat 
transport by advective magma percolation. 

4.4. Viscoelastic Results 

Figure S] compares the tidal work response W(T, ui) 
of the Maxwell, SAS and Burgers models in maps of 
temperature-period space. The figure can be read in hor- 
izontal slices, representing the temperature evolution of 
a planet in a fixed orbit. Or it may be read vertically, 
representing a planet migrating to shorter periods, en- 
countering higher tidal heating as it moves, then shifting 
to higher temperatures in response. The peaks seen in 
each map are a kind of planetary resonance, akin to the 
Debye peak in frequency space. In this case they occur as 
changing mantle temperature in effect tunes the planet's 
response frequency to match its orbital period. Since 
tuning is linear in frequency but exponential in viscosity, 
the peaks appear at roughly constant temperatures. 

Tidal output spans a wide range of magnitudes due to 
the large exponents in equation[6j however Figure[4]is lin- 
ear and clipped to a maximum of 400 TW. This empha- 
sizes the geologic impact of tides, as values much below 
a ~40 TW radiogenic background are of little concern, 
while extreme values above a few hundred TW will move 
the system rapidly in temperature towards global partial 
melting. Note that regions of both negligible (white) and 
extreme tidal heating (black) are quite broad, while the 
geologically moderate range from 10 TW to 400 TW is 
comparatively narrow, suggesting that supcrtidal plan- 
ets only rarely result in enduring Earthlikc features such 
as plate tectonics. In section [5] we show how moderate 



cases may still be favored by equilibrium considerations. 

This figure is drawn for a high e=0.1, where significant 
heating extends out to 40-70 day orbits. Higher planet 
and host masses will extend the viscoelastic range of tidal 
relevance, while lower eccentricities pull the range of tidal 
relevance back within 30 days, as in section [3] 

Unlike the fixed Q approach, Figure 2] also shows that 
declaring a tidal output for a given planet knowing its 
period and eccentricity is not possible without knowing 
its internal mantle state. The range of uncertainly in 
heat due to temperature often spans from negligible to 
extreme. 

On one hand the Maxwell result is similar enough to 
the SAS and Burgers results that it may be sufficient 
to broadly categorize a planet. On the other hand, re- 
sults at a given phase space location can differ by several 
orders of magnitude between models. In Figure |4] full 
rounding of the SAS response peak is expressed since it 
is shifted fully to the cool side of T so i. This is due to 
the SAS model's use of the creep compliance term SJ, 
which is set to 20% below the elastic com pliance J, fol- 
lowing data in Smith ik Carpente r (1987). Still, the run- 
up alone to the Maxwell peak here extends tidal influence 
out further than the SAS model, again due to using J vs. 
5J . Placing a Maxwell peak within our interest range can 
be parametrically difficult without significantly lowering 

Vo- 

The Burgers map is significantly different, with two 
response peaks in the temperature axis. The two peaks 
are well separated with rjg and i]a set to be a factor 
50 apart (e.g. grain boundaries relax 50 x faster than 
voids diffuse). This factor for a bulk planet is uncer- 
tain, w ith laboratory and model values i n the range 2 
to 100 (|Sabadini et al.lfl987l : lPeltieiifl98^ . Smaller val- 
ues p l ace th e response peaks nearer together. Following 
iPostj (|1977h we let E* A =E* B , however small changes can 
significantly alter the peak spacing. 

More advanced models than the Burgers body exist, 
including ones with discrete Voigt-Kelvin subunits in se- 
ries, and the Andradc model, where viscosity and shear 
arc treated as continuum functions in frequency. Re- 
sponse peaks in real materials can be hard to isolate 
(|Cooperll2002fl and real planets will have a range of com- 
positions and grain sizes, blurring an ideal response func- 
tion. In such cases the importance of the Burgers model 
is less the gap between the peaks and more the overall 
response broadening. 

Existing dissipation measurements from Earth's Chan- 
dler wobble, tides and fr ee oscillations, summarized in 
iKarato fc Spetzlerl (|1990D show Q values which change 
gradually with frequency, and are modeled well by em- 
pirical p ower l aws of the form Q ~ co aQ (a p = 0.1 to 0.4) 
(|Karatol [20071 : lEfroimskv fc Lainevl l2007h . Thus, while 
a given Maxwell parameterization may match observa- 
tions well at a 1 year period, it would then typically 
underpredict dissipation at 1-10 days by more than an 
order of magnitude. This is due to the linearity of u> in 
the Maxwell approach, leaving it far more sensitive to 
frequency changes than w 01 or ui 0A . While similar dis- 
crepancies may be expected in other mechanism based 
models, the response broadening in the Burgers model 
does help in alleviating to sensitivity. In general, response 
broadening may be achieved by adding Maxwell elements 
to a model to represent new dissipation modes. For 
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Fig. 4. — Maps showing the magnitude of tidal heating (in TW) as a function of mantle temperature and orbital period for homogeneous 
planets of IMe and e=0.1 around a solar mass host. Viscosity setpoint is 1 X 10 22 Pa-s at 1000K. The solidus and liquidus are set at 1600 
and 2000K respectively. Tidal heat diminishes rapidly above the solidus and is absent above the 1800K breakdown temperature. Broad 
regions of negligible and extreme heating are separated by narrow regions of geologically moderate heat. Gradients are clipped at 400 TW. 



example lO'Connell fc Budianskvl (| 1977ft show that fluid 
flow between cracks and pore spaces of varied size can 
lead to roughly constant Q values over several decades 
in frequency. An Andrade model, either based directly 
on whole Earth data, or based on a continuum model of 
springs and dashpots to model the broad range of compo- 
sitions in the Earth, may overcome many of these issues. 

For this reason, the fixed Q approach of section [3] is 
valuable in worlds that arc particularly Earthlikc. The 
value of derived Maxwell, SAS, and Burgers models 
comes in extrapolating to new bodies, such as super- 
Earths, Hot Earths, Ice Earths, and Earths with var- 
ied compositions, varied partial melt fractions, extensive 
devolatilization, or static lids. 

Many variations of Figure [4] exist. In particular, alter- 
ing the defining viscosity from 1 x 10 22 to 1 x 10 24 Pa-s 
at 1000K shifts the response peaks to above the solidus, 
where they appear sharply compressed along the abscissa 
prior a falloff at T^kdwn ■ Figure [5] helps explain this be- 
havior. 

Figure [5] compares the models as maps in shear- 
viscosity phase space. Overlain on these maps arc tra- 
jectories as a planet evolves in temperature. The up- 
per branch represents the planet prior to melting, when 
shear modulus is constant. Evolving to the left is akin to 
heating with decreasing viscosity. Melting begins at the 
elbow, and the lower partial melt branch represents the 
fall of shear modulus and viscosity together. The slope of 
this branch is only moderately constrained by laboratory 
experiments. 

Beyond the solidus, a given W(T,cu) function will ci- 
ther show a sharply compressed response peak, or a 
dropoff with no peak. These maps help explain why. 
If the downward trajectory crosses a tidal response ridge 
or island, a sharp post-solidus peak will occur. If the 
downward trajectory slips off the side of an island, or 
passes between the saddle of the two Burgers islands, no 
post-solidus tidal peak occurs. 

Both Figure 2] and [5] support the notion that a broad 
range of models and parameters result in significant 
phase space volumes where tidal heating is extreme, co- 
existing with broad negligible regions. We have explored 
variations of the Burgers parameters t]a and t)b to verify 
they move and stretch the dual features of the Burgers 
plots, but preserve the function's basic form. Variations 



in the creep activation energy E* (A and B for Burg- 
ers) has a strong impact on final magnitudes, but again 
smoothly deforms the functions. 

Because of the overall uncertainty in both parameters 
and models, we find it misleading to report a specific vis- 
coelastic tidal output for a given planet without detailed 
analysis. Of more importance is the broad range of out- 
puts possible. However, due to feedback with convection, 
a planet is not equally likely to reside at all possible tem- 
peratures, and instead will seek out specific equilibria. 

5. TIDAL EQUILIBRIA 

Following the method of IMoord (|2003allbh for Io, we 
can model the thermal behavior of a supertidal body in 
time by simultaneously plotting tidal work as a function 
of temperature W(T) atop convective output E conv (T). 
Equilibria points occur whenever tidal input equals con- 
vective output. 

We use the method of parameterized convection 
(jO'Connell fc H agcr 1980) to iteratively solve for cooling 
as a function of mantle temperature. Convective vigor 
determines the thickness of a conducting boundary layer 
S by the formula: 



d 
2a~2 



Ra 
Ra c 



(22) 



a%: Flow geometry constant ~1 

d: Mantle thickness ^3000km 

Ra c : Critical Raylcigh Number (free top)'- 



4100 



where the exponent —1/4 is true for bodies with in- 
ternal heat generation (—1/3 for bodies without). The 
Rayleigh number is found from the material properties of 
the uniform interior and heat flow through the top qbl'- 



Ra 



a g p d 4 



f]{T) K kthe 



(23) 



g: Gravitational acceleration ~9.81 m/s 2 
a: Thermal expansivity, ^lxlO~ 4 (|Kaulalll968f ) 
p: Average Density, ^5000 kg/m 3 
At: Thermal diffusivity = k t herm/ pC p 
ktherm'- Thermal conductivity, ~2 W/mK 
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Fig. 5. — Maps showing the magnitude of tidal heating as a function of shear modulus and viscosity for homogeneous planets of IMe, 
e=0.1, and a period of 10 days. Melting trajectories are overlain, and breakdown is marked with a dot. Changing the orbital frequency 
shifts the underlying map relative to given fixed melting path. The generally unsuitable Voigt-Kelvin model is included for comparison. 
Note how in the Burgers model a cooling planet can cross a radically different history of islands than in the Maxwell or SAS models. 



Because of the viscous term rj(T), the Raylcigh Num- 
ber responds rapidly to temperature. We neglect the 
weaker temperature dependencies in p and a. Conduc- 
tion through the top boundary is computed based on 
an average thermal conductivity htherm an d the thermal 
gradient. 



Qbl 
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An initial guess is made for boundary layer thickness, 
then the model is iteratively solved for qbl(T). We as- 
sume the planet is in equilibrium with its star and has 
a thin atmosphere, allowing all heat that moves through 
the boundary layer to escape into space. 

Knowing the shape of W(T) (plus a small radionu- 
clide background) and E conv (T) essentially characterizes 
the behavior range of a supertidal planet in time. The 
crossing points of the two functions represent equilibria 
points. These may either be stable or unstable based on 
the relative values of dW/dT and dE conv /dT. Figure [6] 
shows an example for mild SAS tides. Figure [7] shows 
a stronger example for a Burgers body. Young planets 
subject to tidal heating early may evolve into this system 
from high temperatures, while planets that significantly 
cooled prior to onset tidal forcing approach equilibria 
from the left. 

When tidal forcing is strong a hot stable equilibrium 
point typically exists near the onset of melting, where 
viscosity falls off sharply. Tides and convection cross 
here, with tides shutting down just as convection sharply 
ramps up. The location of this equilibrium is robust in 
cases when it occurs prior to onset melting. In cases 
where it occurs after melting onset, due to our homoge- 
neous mantle assumption, only the existence of the equi- 
librium is robust. The predicted degree of partial melt \ 
is at best suggestive of how far a real mantle will melt. 
Once at a stable equilibrium, a planet will remain there 



at constant bulk temperature as long as the orbit allows. 

If tides are weak, W(T) may never intersect E conv (T), 
meaning no equilibria exist, and secular convective cool- 
ing will proceed uninterrupted, although perhaps at var- 
ied rates. In special cases (see Figure [6]) heating may 
nearly be enough to halt secular cooling, slowing cooling 
down to a near standstill for 10-100 million years, before 
more rapid cooling resumes. 

When a planet evolves across a response peak in either 
direction it experiences a wave of internal heating. Inter- 
nal temperature changes are in effect tuning the planet's 
natural frequency to match a constant orbital frequency. 
In the Maxwell and SAS cases, a single episode of sudden 
peak heating can occur (regardless of whether the heat 
is more or less than convection). In the Burgers case 
two episodes are possible. Changes in orbital forcing, or 
thermal-orbital feedback could cause planets to evolve 
though these events multiple times. 

The time it takes a planet to reach equilibrium depends 
mainly on initial conditions. Our simulations show typi- 
cal tidal response peaks are crossed rapidly, on the order 
of 10-50 million years. This will be manifested in the 
planetary history as a sudden episode of extreme heat- 
ing, possibly recorded on the planet's surface, followed 
typically by more moderate equilibrium heat rates. 

We looked for cases where the peak in W(T) could lead 
to cyclic overshoot events but found the system dynami- 
cally ovcrdamped, with cyclic, quasiperiodic, and chaotic 
solutions prevented by a planet's high thermal inertia 
and long heat transport timcscalc. Single overshoots do 
occur, in particular after heating across a strong reso- 
nance peak when the hot stable equilibrium is well below 
the solidus. 

While equilibrium considerations help limit the range 
of likely states, uncertainty remains. Extreme equilib- 
rium tidal heating can occur, if by parametric chance the 
convection curve crosses exactly at the response peak, 
creating a stable equilibrium very near the maxima in 
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Fig. 6. — Behavior Example 1. SAS Model. Nearly halted secular cooling. Panel (a): Tidal work profile (solid line) vs. convective output 
profile (dashed line) (e=0.044, lMg, 20 day period, material properties from Tablc[5] hot start, negligible radionuclides, 0.5Mg o ( primary). 
Panels (b) and (c): Bulk mantle temperature and tidal heat evolution in time. Eccentricity is set for peak tidal heating to almost match 
convective cooling. Note the long timescale and long plateau as the system slowly evolves over the response peak. This forms one limit of 
slow secular cooling. Higher tidal forcing leads to tidal-convective stabilization and a trivial equilibrated time history. Lower tidal forcing 
leads to a weaker single episode of heating. Cold initial conditions will cool without peak crossing. 
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Fig. 7. — Behavior Example 2. Burgers Model. Panel (a): Tidal and convective profiles. Hot stable equilibrium marked with solid circle. 
(e=0.05, IMe, 12 day period, material properties from Table[2] cool start, lMg a i primary). Panels (b) and (c): Dual peak crossing: this 
strong forcing case heats the planet through both Burgers body resonance peaks, such that two unequal sudden warming episodes occur in 
time. Note the short timescale. Peak heating is ~12200 TW while equilibrium occurs at ~580TW and 1646K (Bulk x=H-5%). Hot start 
planets would cool to the hot stable equilibrium without encountering the peaks of this scenario. 

W(T). With favorable parameter choices, such extreme 
tidal equilibrium solutions can still rise as high as mil- 
lions of terawatts, as in the fixed Q approach. Therefore, 
as with the earlier extreme solutions, we need a better 
understanding of localized partial melting and advective 
cooling to know if such high outputs occur in real planets. 

We can generalize behaviors further to incorporate 
changes in orbital forcing by constructing a bifurcation 
diagram as shown in Figure [8] Movement of equilibria 
points arc shown here as functions of eccentricity. Sim- 
ilar diagrams can be created for variations of the semi- 
major axis. These diagrams are in effect formed by tak- 
ing a tidal heating curve W(T) and moving it without 
deforming its shape vertically past the convection curve 

Econv ■ 

Bifurcations occur as the peaks in W(T) make first 
contact with the convection curve (see Figured). At low 

eccentricities, W(T) lies entirely below E conv (T), so tidal 
heating is weak, never competing with cooling, and tidal 
equilibria are impossible. Bifurcation points represent 
threshold eccentricities where tidal equilibria come into 
existence. 

A stable branch, or locus of stable equilibria points, is 
shown as the solid line. The hot stable branch typically 



starts near or just after melting and shifts gradually to 
higher degrees of partial melt. An unstable branch forms 
on the cold side, shown as the dashed line. Below the bi- 
furcation point, eccentricity is too weak and only secular 
cooling is possible. 

The region above the two branches represents heating, 
while the region below represents cooling. Trajectories of 
planets in the heating region evolve towards the hot sta- 
ble branch, becoming trapped there at equilibrium. Tra- 
jectories evolve away from the unstable branch. Rates 
of evolution are based on W(T)-E conv (T). Time varying 
radionuclides are not included Figure [8l but will gener- 
ally add a quasi-stable state at lower temperatures. Or- 
bital migration or eccentricity pumping can move trajec- 
tories vertically from the cooling regime into the heating 
regime or vice versa. Eccentricities in resonances oscil- 
late around a fixed value, which would be represented 
by a sinusoidal trajectory in this diagram, and has the 
potential to shift marginal planets between the heating 
and cooling regimes in a complex manner. 

What these diagrams best demonstrate is that plan- 
ets will either tend to evolve quickly to the hot stable 
branch, or will never face tidal heating that can compete 
with convection and will secularly cool, although perhaps 
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Fig. 8. — Bifurcation diagrams of tidal equilibria for a IMe planet in a 16 day orbit, comparing the Maxwell, SAS, and Burgers models, 
each for two rheologies. Loci of stable equilibria: solid lines. Loci of unstable equilibria: dashed lines. Regions above the equilibria 
branches represent planetary heating (tidal heat exceeds convection). Regions below the branches represent cooling. Cooling from a high 
melt fraction will come to a halt on the stable branch if tidal forcing (eccentricity) is high enough. Panel (b) shows a hypothetical trajectory 
of migrating into a high eccentricity state. Panels (a), (b), (c): baseline rheology. Panels (d), (e), (f): dry or devolatilized rheology with 
higher defining viscosity, where response peaks shift above the solidus. For some Burgers cases the low temperature stable branch may 
dominate, particularly at longer periods. 



at a slower rate than without tides. This supports our 
discussion in section [2^4l that maximum dissipation states 
tend be brief, unless equilibrium occurs at a W(T) peak. 

Figure [5J:,f shows the same information for a Burg- 
ers body with two response peaks, thus two bifurcation 
points, two stable branches, and two unstable branches. 
More complex planetary histories may occur. In particu- 
lar, planets may become trapped at a colder tidal equilib- 
rium associated with the grain boundary slip mechanism. 
This cold stable branch is even more pronounced at peri- 
ods >20 days and can be the first stable state reached by 
many planets. However, as inhomogencous mantles may 
blur distinct peaks, our Burgers results are best viewed 
as a demonstration of the increase in behavioral com- 
plexity that occurs when additional response frequencies 
are taken into account. 

6. DISCUSSION 

This work highlights the question of what will be the 
ultimate shutdown mechanism for an extreme tidal ter- 
restrial planet. Both the fixed Q method and the gener- 
ally more conservative viscoelastic methods predict that 
in some circumstances tidal heating can reach millions of 
terawatts within a planet modeled as homogeneous. Our 
models of tidal-convcctivc equilibria are very effective in 
exploring planetary behaviors prior to equilibration, but 
only coarsely resolve actual equilibrium heat rates due to 
the assumption of homogeneity. 

6.1. Inhomogeneous Melting 

Onset partial melting can begin in an inhomogeneous 
planet at much lower heat rates. To rou ghly determine 
the lo cation of melt initiation, we follow TValencia et alj 
(|2006h to calculate the temperature profile with depth, 
or geotherm T(z), of a tidally heated exoplanet. We esti- 
mate the solidus as a function of pressure via the Simon 
Law us i ng da ta for forsterite, enstatite, and iron from 
iPoirierl (|2000f ) (neglecting high pressure mineral phases 



and eutectic mixtures) . While some models of Earth sug- 
gest partial melt already exists at the lithosphere's base, 
we find the geotherm can begin to intersect the solidus 
with as little as 30-40 TW of added tidal heat, an ap- 
proximate doubling of Earth's current output. Where 
sufficient local partial melting occurs, tidal friction will 
decrease while continuing in better tuned viscous regions. 
This suggests how a planet may have difficulty generating 
the millions of TW solutions found earlier in this paper. 

We also find core temperature is a nearly linear func- 
tion of tidal input, primarily because of the strong lin- 
ear dependence of the conductive geotherm through the 
lithosphcrc on total heat flow. We assume no tidal heat 
is deposited in the core itself, however small amounts of 
tidal heat (<10TW) can shift the geotherm such that 
the entire core becomes liquid (based on the shallow 
slope of the Simon Law solidus for pure Iron). Thus 
even weak tidal exoplanets may have no inner cores, dis- 
rupting magnetic dynamo activity, just as with younger 
exoplanets prior to core crystallization. 

Models were tested with an upper-lower mantle ther- 
mal separation at 660km depth, as well as asthenosphcric 
only tidal input. Onset melting results were largely the 
same, except for a weaker dependence of the core temper- 
ature on tides, since vigorous upper mantle activity led 
to thinner conductive lithospheres. Varying the planet's 
mass, we find above ~2.QMe the mantle adiabat may 
curve sufficiently for onset melting to also occur at the 
core-mantle boundary (using a Birch-Murgnahan equa- 
tion of state). 

As tides increase, we expect partial melt regions to 
grow in volume, robbing the mantle of material that is 
properly viscoclastically tuned, and bringing about the 
process of tidal shutdown. This requires us to abandon 
a singlc-gcotherm model, since the role and volume frac- 
tion of both rising and falling convection plumes becomes 
key. In the rising plumes, the geotherm may remain 
above the solidus in much of the upper mantle. In the 



Viscoelastic Tidal Exoplanets 



15 



falling plumes, the geotherm may remain entirely sub- 
solidus. Tidal heat shutdown will be governed by what 
volume fraction of the mantle exists within higher partial 
melting zones, and what volume remains best tuned to 
the orbital period. For certain rhcologies, ideal tuning 
may occur in x=l-2% regions, while full detuning occurs 
at only a few percent greater melt. 

For the most vigorous exoworlds, the assumption of 
convective heat transfer gives way t o a regime dom- 
inate d by advective heat transport ([Spiegelman et aTl 
120011 ). Percolation of magma to the surface is vastly 
more efficient than either bulk convection or conduction, 
and can move heat rapidly through the bottleneck of a 
lithosp here. This is akin to the Heat Pipe mod el sug- 
gested (|Moorell200lHMonnereau fe Dubuffetll2002T ) to ex- 
plain Io's apparently strong lithosphere and tall moun- 
tains. Introducing advective cooling in Figure Ek or [TK 
is equivalent to the E conv (T) curve tipping even more 
sharply to vertical as soon as melts begin to mobilize. 
This alone does nothing to set an ultimate limit for tidal 
heat production. Percolation in fact allows much greater 
tidal heat rates, since the mantle can stay predominantly 
rigid, meanwhile running an efficient fluid cooling net- 
work, suggesting that tens of millions of terawatt tidal 
solutions remain geologically possible. 

The true internal extremes of tidal heat production 
will be governed by models that include grain size and 
porosity estimates. The amount of cracks and punctures 
in a lithosphere and the ability of flow networks to en- 
dure in a shifting mantle will become vital. An imperme- 
able mantle leads to high melt storage, and a reduction 
in tidal behavior. An easily permeable mantle leads to 
an efficient melt network with sparse storage and broad 
tuned- viscosity zones that allow significant tidal coupling 
to continue. Vigorous tidal processing may lead to sig- 
nificant global chemical layering between refractory and 
low melting temperature species. 

6.2. Magma Oceans 

Hot Earth planets may have insolation supported 
magma oceans where basal friction due to tidal slosh 
plays an important role. However mantle tidal heating is 
unlikely to initiate exposed surface magma oceans from 
the top down. Fresh lakes of Haw aiian lava solidify t o 
a depth of ~8 meters in one year (IWright et al.lll976f). 
For a heat of fusion near 506 kJ/kg (jKoiitani fe Akaogil 
Il995t ). tidal heating can advectively resurface a IRe 
planet at a rate of ^0.0024 cm/yr/TW (given Earthlikc 
background heating). At this rate, 330,000 TW of tidal 
heat are needed to resurface the entire planet at 8 m/yr, 
to begin sustained burial of unsolidificd material (Assum- 
ing all lava flows start at T so i and planetary heat is lost 
solely through magma volcanism). Full magma ocean 
build-up via outpourings is likely to require half a mil- 
lion TW or more. So while planet-wide resurfacing may 
occur, it is difficult for tidal heating to build up a surface 
magma ocean with no assistance from insolation. Alter- 
natively, sub-lithosphcric melting and subsequent thin- 
ning may produce insulated near-surface magma oceans 
or crystalline slush layers at lower tidal rates. 

Individual volcanic vents may produce lava flows 
greater than 8m/yr and produce significant localized 
magma lakes. For example, only 370 TW of tidal heat 



can resurface 1% of Earth in lm deep lava flows per year. 
Thin-layer global resurfacing as on Io is unlikely for vis- 
cous lavas. This supports the notion of searching for 
small radiantly cooled hotspots on supertidal exoplan- 
ets. 

Self gravity in W(T, ui) raises the possibility of a super- 
tidal magma ocean planet. Direct fluid tidal work is very 
small, due only to the low fluid viscosity. Self-gravity 
allows for greater peak outputs through resonance by 
adding an effective ~ 2 x 10 11 Pa spring. For a fluid 
planet: 



Im(k 2 ) 
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Response timescales for a fluid Earth are on the or- 
der of seconds even for the thickest conceivable magmas 
{VHq = 0-01 to 1000 Pa-s depending strongly on composi- 
tion, pressure, and volatile content (|McBirnev fc Murasel 
119841 )). Still, the climb toward resonance at very short 
periods does occur in the fluid regime just beyond 
Tbrkdwn- Fluid tidal heat may play a role for terrestrial 
planets below 1 day orbits whose mantles never crystal- 
lized, or whose insolation induced magma oceans have 
grown to great depths. 

Following Murray and Dcrmott, surface deformations 
of l-10km are possible on extreme tidal bodies in 1-2 
day orbits (deformation is only of the order 300m for 
Io, and 20m for Europa). This suggests massive faulting 
and cracking if a lithosphere were to remain brittle, or 
powerful tidal currents if material is melted. Note e ma x 
is a function of the orbit only, not on Q or Eudai- The 
structural geology of extreme periodic fault excitation is 
particularly interesting in light of the tiger-stripe stress 
concentrations on Enceladus. 

6.3. Habitable Zone Modifications 

Finally, we briefly consider tides and habitability. The 
conventional definition of habitable zone (HZ) uses lumi- 
nosity to set a range of semi-major axes for stable liquid 
surface water ((Easting et al. 1993:!) While chemistry and 
outgassing rates do matter for true habitability, we here 
consider temperature only. 

Tidal heating primarily provides additi onal habitabi l- 
ity opportunities for extrasolar moons (|ScharJ [2001 . 
Tides on larger eccentric moons may not only generate 
liquid subsurface oceans as with Europa, but also liquid 
surface water, all at arbitrary distances from a host star. 
To demonstrate, in Figure [5k., we explore a range of or- 
bital parameters for hypothetical 1 and 7 Me moons in 
orbit around a lMj host, using equation [6] to estimate 
tidally driven blackbody surface temperatures. For such 
large moons, 273K<T stlr /<373K is possible for many 
reasonable values of period and eccentricity given zero 
contribution from insolation. Large moons of this type 
may be rare, however we are being conservative looking 
only at surface temperatures. Subsurface oceans can be 
sustained at more relaxed masses and eccentricities. Host 
masses above lMj also enhance tides. Thus provided an 
enduring eccentricity source exists, tides may support is- 
land habitable zones around nonluminous p rimaries such 
as ou ter planets, pulsar s, ultracool dwarfs (|Martm et al.1 
Il999f ). isolate planets ijOsorio et al.l l2000h . and ejected 
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Fig. 9. — Habitable Zone Modification. Panel (a): Habitable 
zones for a 1 and 7Mg moon orbiting a nonluminous lMj host. 
Blackbody temperatures for liquid surface water are achieved in 
orbits similar to the Galilean system with elevated eccentricities. 
Panel (b): Habitable zones break up into a spectrum by mass at a 
M3V dwarf of 0.001L SO ;. Panel (c): Habitable zones at a brighter 
dwarf of 0.01L SO ;. Here the IMg HZ is nearly unaltered by tides, 
but the super-Earth is large enough that its inner HZ edge is pref- 
erentially reduced. All panels use SAS tides with a near-melting 
isoviscous mantle at 1 X 10 17 Pa-s, akin to Q=50 ^2=0.3, and an 
Albedo of 0.3. Panels (d), (e), (f): Habitable Zone width reduc- 
tions as a function of eccentricity for the cases in panels (a), (b), 
(c). Widths generally contract until tidal heat exceeds insolation 
heat. 

planets t hat either retain or accum ulate large resonant 
satellites (jDebes fc SigurdssonlfeOOTf ). 

A second aspect of tides and habitability is how tidal 
heat modifies existing habitable zones. Tides can break a 
HZ into a spectrum by mass, as in Figure [SJx In general 
tides will also reduce HZ widths, due to the a~ 6 term in 
equation [6l A given increment in eccentricity will lead 
to a greater change in surface temperature at the inner 
edge of a HZ than at the outer edge. Width reduction 
occurs rapidly at low eccentricities, as in Figure^). Some 
width recovery occurs once tides become the dominant 
heat source, but this effect will be rare due to the high 
eccentricities required. Our calculations do not include 
the effect of average eccentric insolation, and thus we do 
not calculate beyond e=0.5. Width reduction is more 
rapid for super-Earths, and can occur preferentially in 
special borderline cases as in Figure OU. 

At stars brighter than ^O.OlLsoi habitable zones will 
be unaltered by tides. One can think of a region in semi- 
major axis space akin to a habitable zone called a tidal 
zone (TZ), a region where tides can influence surface tem- 
perature, typically very close to a star. For G and K 
stars, the habitable zone is much further out than the 
tidal zone, and tides never alter habitability. For M3V 
dwarf stars and dimmer, the HZ and TZ begin to over- 



lap. At this transition, preferential HZ reduction occurs 
at the inner edge (for 1 vs. 7 Me). At lower luminosi- 
ties, a HZ is primarily shifted outwards with less width 
reduction. While dwarf stars such as GJ876 (0.0124L so z) 
have not been traditional targets for life searches, inter- 
est has risen since moons of gas giants can overcome the 
problem of spin-synchronization, and ice shells or strong 
magnetic fields can improve radiation shielding. In all 
cases we have used an SAS planet with a near-solidus 
isoviscous mantle of 1 x 10 17 Pa-s, what may be consid- 
ered ideal viscoelastic tuning. These results are similar to 
using Q=50 and &2=0. 3 and thus simi lar to the measured 
Earth response values summarized in iKarato &: Spetzlerl 
(|1990f h Cooler or dcvolatilizcd planets may be too stiff 
for tides to couple this effectively. 

7. CONCLUSIONS 

In this paper we have shown how a range models pro- 
duce extreme tidal heating in short period terrestrial ex- 
oplancts. The existence of broad regions of extreme tidal 
solutions lying alongside negligible solutions is robust to 
parameter uncertainty. However this dual nature makes 
it difficult to specify a given planetary heat output based 
on tidal forcing strength alone, without knowledge of the 
interior. Broadly we find tidal heating in excess of ra- 
dionuclide heating occurs below approximately 10-30 day 
orbital periods. Tidal heating in excess of insolation oc- 
curs only below ^2 day orbits, and may not be realized 
due to onset melting and tidal decoupling. Tidal heating 
around M dwarf stars has the best chance to match or 
exceed insolation. 

Due to the dominance of insolation at short periods, 
tidal heating has a negligible impact on the blackbody 
surface temperature of planets hosted by K class stars 
and brighter. For such worlds, tidal relevance is pre- 
dominantly internal. When insolation is strong, moder- 
ate temperature elevations due to extreme tides of 1-5 
degrees may easily be masked by a small percent un- 
certainty in albedo. Observational confirmation of ex- 
treme tides on a candidate planet suggested by eccen- 
tricity should therefore focus primarily on chemical and 
hotspot signatures. 

Our viscoelastic parametric studies show broad regions 
in phase space where short period exoplanet tidal heat- 
ing is both geologically negligible (<1 TW), and extreme 
(>400TW), with a narrow band of moderate tidal cases 
(1-400 TW) separating the regions. The Burgers body 
and SAS models can reveal otherwise hidden features in 
tidal evolution not manifested by the simpler Maxwell 
model, including extra episodes of warming and extra 
equilibrium points. Use of the creep compliance term 6J 
in place of the bulk compliance J can lead to significant 
variations in results. The distinct response peaks that 
occur in a Burgers model with well defined parameters 
may however be blurred in a real planet with a range of 
heterogeneities. 

If a planet can ever couple significantly with tides, it 
will either secularly cool into a stable equilibrium with 
convection, or undergo 1-2 rapid warming episodes be- 
fore stabilizing. Equilibrium heat outputs are sensitive 
to model choice, can be negligible or extreme, and can 
occur both above or below the solidus. Tidal heating 
typically reaches equilibrium with convection in a few 
million years. Stable equilibria can be shifted, created, 
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or destroyed by changes in eccentricity. 

Melt advection is the key to determining the limits 
of terrestrial tidal heat production. True tidal equilibria 
must be derived from a balance between partial melt pro- 
duction regions and melt percolation rates, invoking sep- 
arate adiabats for rising and falling convective plumes. 
Tides are unlikely to generate surface magma oceans on 
their own, but due to high insolation at short periods, 
they may help control the depth of preexisting magma 
oceans. 

Extreme tides, if present, will mainly alter the hab- 
itable zones for short period planets of lower luminos- 
ity M dwarf stars. Habitable zones may be shifted out- 
wards, preferentially reduced from their inner edges, and 
partly or completely split into a spectrum by mass. Tides 



at brighter stars will only influence uninhabitable hot 
worlds, or may alter cold resonant moon systems at ar- 
bitrary distances from the host star. 
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